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Abstract 



We investigate whether present data on helioseismology and solar neutrino fluxes may constrain 
WIMP-matter interactions in the range of WIMP parameters under current exploration in WIMP 
searches. We find that, for a WIMP mass of 30 GeV, once the effect of the presence of WIMPs in 
the Sun's interior is maximized, the squared isothermal sound speed is modified, with respect to 
the standard solar model, by at most 0.4% at the Sun's center. The maximal effect on the 8 B solar 
neutrino flux is a reduction of 4.5%. Larger masses lead to smaller effects. These results imply 
that present sensitivities in the measurements of solar properties, though greatly improved in recent 
years, do not provide information or constraints on WIMP properties of relevance for dark matter. 
Furthermore, we show that, when current bounds from direct WIMP searches are taken into 
account, the effect induced by WIMPs with dominant coherent interactions are drastically reduced 
as compared to the values quoted above. The case of neutralinos in the minimal supersymmetric 
standard model is also discussed. 
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I. INTRODUCTION 



A host of independent astronomical observations point to the existence in our Universe of 
a total amount of matter in the range 0.2 < tt m < 0.4 (or equivalents, 0.05 < VL m h 2 < 0.3, 
where h is the Hubble constant in units of 100 km • s _1 - Mpc -1 ), well beyond the amount 
of visible matter fl vis ~ 0.003. Since the primordial nucleosynthesis tells us that baryons 
cannot contribute for more than about 5% 0, most of the dark matter must be non-baryonic. 
The evolution theory of the primordial density fluctuations into the present cosmological 
structures indicates that most of the dark matter must be comprised of cold particles, i.e. 
of particles that decoupled from the primordial plasma when nonrelativistic. A particle with 
the suitable properties for being a significant cold dark matter relic is generically defined as 
a WIMP (weakly interacting massive particle). 

A variety of physical realizations for a WIMP are offered by various extensions of the 
Standard Model ||, the neutralino being one of the most appealing candidates. In the 
present paper, many considerations are developed in terms of generic WIMPs. We consider 
the neutralino, whenever we wish to narrow down to a specific candidate. 

WIMPs are very actively searched for by means of various experimental strategies. Direct 
searches rely on the measurements of the signal that a nucleus of an appropriate detector 
would generate, when hit by a WIMP Indirect searches are based on measurements of 
the signals due to WIMP pair annihilations in the galactic halo or inside celestial bodies || . 

The event rates of WIMP direct search experiments are proportional to the product of 
the local WIMP density p x in the Galaxy times the WIMP-nucleus cross-section. In what 
follows p x will be expressed as p x — £ • pi, where pi is the local total dark matter density and 
£ (£ < 1) is a scaling parameter which accounts for the actual fraction of local dark matter 
to be ascribed to the candidate x- 

WIMP-matter interactions are conveniently classified in terms of coherent cross-sections 
and (nuclear) spin-dependent ones. In the first case, the WIMP-nucleus cross-section is 
given in terms of the WIMP-proton and WIMP-neutron cross-sections. To simplify the 
formalism, in the following we further assume that, in the coherent case, WIMPs interact 
with equal strength with protons and neutrons (for instance, this may safely be assumed for 
neutralinos, while it is not the case for neutrinos), and then we express a generic WIMP- 
nucleus coherent cross-sections in terms of a single WIMP-nucleon cross-section a c . For 
the spin-dependent case, the derivation of a WIMP-nucleon cross-section from the WIMP- 
nuclear one depends on the nature of the WIMP and on specific nuclear properties || . 

In the case of WIMPs with coherent interactions with matter, the range of WIMP pa- 
rameters under current exploration in WIMP direct searches is conveniently expressed in 
terms of the quantity £cr c : 

10~ 43 cm 2 < £a c < 6 • 10" 41 cm 2 , (1) 
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for WIMP masses m x in the range: 

30 GeV < m x < 270 GeV. (2) 

In the derivation of the ranges in Eqs. one has taken into account a variety of 

WIMP distribution functions in the galactic halo and uncertainties in the determination of 
the relevant astrophysical quantities || [7J. 

For the reasons mentioned above, in the case of spin-dependent interactions no model- 
independent sensitivity range may be derived. However, to give an indication, we recall 
that, in the case of a relic particle which interacts with matter mainly by spin- dependent 
interactions mediated by Z-boson exchange, current WIMP direct experiments are sensitive 
to values of ^a s d of about 10" 36 — 10~ 37 cm 2 ||, where a S( i is the spin-dependent WIMP- 
proton cross-section. 

One of the WIMP direct search experiments, the DAMA/NaI(Tl) experiment ||, has 
observed an annual modulation effect (at a 4<r C.L.) with all the features expected for a 
WIMP signal |J. When interpreted as due to a WIMP with coherent interactions, the 
DAMA effect provides a 3a region in the plane m x — £<r c embedded in the range of Eqs.(jl|- 
0). In Ref. |IJ it was proved that this annual modulation effect is compatible with an 



interpretation in terms of relic neutralinos (see also Ref. [|TT| ) . 



If WIMPs populate our Galaxy, they can be captured by the Sun and accumulate in 
its core. For suitable values of the WIMP parameters, this could have relevant effects on 
the Sun. In fact, WIMPs could provide an effective mechanism for energy transport in the 
Sun, producing an isothermal core and reducing substantially the Sun central temperature 
T c . Following this idea, some time ago a special class of WIMPs named cosmions, with 
masses of a few GeV and scattering cross sections on nucleons of the order of 10 -36 - 10~ 34 
cm 2 , was studied in detail as a way of solving simultaneously the solar neutrino puzzle and 
the dark matter problem (see e.g. Refs. |12|, [13|, [14], [15], [16], [T7|])- The cosmion hypothesis 
was progressively abandoned when it became clear that the solar neutrino puzzle cannot be 
accounted for by simply reducing T c . 

In the last twenty years, our observational knowledge of the solar interior has progressed 
enormously. By means of helioseismic data it has become possible to derive the sound speed 
with an accuracy of about one part per-thousand over most of the radial profile and of 



about one percent in the innermost part [fL8| . By the same method, it has been possible to 
deduce important properties of the convective envelope. The photospheric helium fraction 
Y p h and the depth of the convective envelope Rb have been determined with an accuracy of 
about one per cent and one per-thousand respectively, following the pioneering papers of 



Refs. [19, 20]. Moreover, the measurement of the neutrino flux from B decay obtained by 



combining |^T], |22|, ^3[ |24| SNO charged current |25|] and Super-Kamiokande []26|] data and 



confirmed by the recent SNO neutral current results [27], has provided a determination of 



the temperature T c near the center of the Sun at the level of about one per cent [28 . 
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All the predictions of the Standard Solar Model (SSM) have been confirmed by these 
accurate tests, so one is naturally led to the question of whether our accurate knowledge 
of the Sun interior can be used to constrain the WIMP parameter space. This particularly 
interesting question was recently raised in Refs. [29| and [30], where it was concluded that 
solar physics can be used to significantly constrain the WIMP parameter space. 

In this paper, we provide an alternative analysis of the problem, obtaining results sub- 
stantially different from those derived in Refs. |29], |30|. 

The plan of the paper is as follows. In Section II we discuss WIMP energy transport 
in the Sun, introducing the relevant physical parameters to be discussed in the subsequent 
sections. In Section III and IV WIMP energy transport is applied to the central solar 
structure, in order to determine the regions of the WIMP parameter space where the WIMP 
energy transport is most efficient. This leads to an optimal choice of WIMP parameters 
which is then used in Section V to calculate the solar sound speed profile and discuss the 
possible impact of helioseismology on WIMP properties. The sensitivity of the available 
data on 8 B solar neutrinos to the presence of WIMPs in the Sun is then discussed in Section 
VI. The neutralino, as a specific realization of WIMP, is discussed in Section VII. Section 
VIII is devoted to our conclusions 



II. WIMP ENERGY TRANSPORT 



In this section, we briefly review some issues relevant to the problem of WIMP energy 
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transport in the Sun addressed in Refs. [|12| |13| , |14| , |15| , [U 

WIMPs are captured by the Sun during the course of its lifetime and are confined by the 
effect of gravity in a central region with a radial number density distribution approximately 
given by pi], R2f (we hereby use natural units: h — c — ks — 1): 



n x (r) = {n 



cxp 



(3) 



where the quantity 



the total number of WIMPs in the Sun and 



(1/7T 3 / 2 ) (N x /rl) is the central WIMP number density, N x is 



3T r 



1/2 



\2irG p c m^ 



(4) 



T c and p c being the Sun central temperature and density, and m x the WIMP mass. For 
T c = 1.57 ■ 10 7 K and p c = 154 g cm" 3 , one finds r x ~ 0.01 (100 GeV /m x ) l l 2 R Q . We note 
that Eq. (^) is obtained under the assumption of constant density and temperature in the 
region of the Sun populated by WIMPs. This approximation is accurate enough for our 
purposes. 
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WIMPs provide a mechanism for energy transport in the solar core whose efficiency 
depends on the Knudsen parameter 

K=!&, (5) 
r x 

where / x (0) is the WIMP mean free path near the solar center and generally: 

= 5>X,(r) (6) 

Here the sum extends to all the nuclear species present in the Sun, is the mass fraction 
of the z-th element, p(r) is the density profile of the Sun and <jj is the WIMP scattering 
cross section with the nucleus of species i. 

For X<1 the WIMP mean free path is much smaller than the dimension of the region 
where the WIMPs are confined and the energy transfer mechanism is conductive. The 
problem of thermal conduction by a dilute gas of massive particles is discussed in Ref. |L2] . 



For our purposes, it is useful to express the WIMP luminosity in the conductive regime by 
using the radiative transport equation: 

167ra rad T(r) 3 r 2 dT 
x ' condlrj 3 K x (r)p(r) dr } 1 ] 

where the WIMP opacity k x is defined as: 



3n x (r) 



K x (r) 4a rad T(r) 3 \ 



' (r) P (r)l x (r)X x . (8) 



In Eqs. ([7|,|]) the quantity a ra d denotes the radiation density constant and X x is a dimen- 
sionless coefficient which depends on the ratio between the WIMP mass and the masses of 
background nuclei 1 . 

For K 3> 1 the WIMP mean free path is much larger than the region in which WIMPs 
are trapped and the energy transfer mechanism is non-local. An analytic approximation to 
treat this case has been developed by Spergel and Press |13j who described WIMPs as an 
isothermal gas at an appropriate temperature T x . In this case, the luminosity £ x ,nl carried 
by WIMPs can be expressed as: 

L x Mr)= f dr' 4nr' 2 p(r') e SP (r';T x ), (9) 
Jo 



The coefficient A x can be derived from the values of the thermal conduction coefficient K(m x /nii) tabulated 
in Ref. JlJ]. It holds: A x = (l/'x)[Ei(lAx,i where is the WIMP mean free path relative to the 

nucleus i. 
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where esp represents the energy transferred to WIMPs per second and per gram of nuclear 
matter and is given by: 

e SP (r;T x ) = 8 J^n x (r) x 

Mr) rn^rn / ' m x T{r) + m 8 T x \ 1/2 
Y 4 ™i K + rrii) 2 \ m x mi J 

The WIMP temperature T x is fixed by requiring that the energy absorbed by WIMPs in 
the inner region (where T(r) > T x ) is equal to that released in the outer region (where 
T(r)<T x ),i.e.: 

L x MRo) = [ R& dr' 4vrr' 2 p(r') e S p(r'; T x ) = . (11) 

J 

Gould and Raffelt | I4]| , by means of a Montecarlo integration of the Boltzmann collision 
equation, showed that the Spergel and Press approximation, Eq. (|10D , overestimates the 
energy transfer typically by a factor of a few. Moreover, by studying simplified stellar 
models, they showed that it is possible to approximate energy transfer in the non-local 
regime by a conductive treatment if one applies a global "luminosity suppression factor". 
We follow their approach by using as a general expression for the WIMP luminosity: 

L x (r) = f(K) L x , cond (r) , (12) 

where 

«^(WTT' (13) 

and K = 0.4 is the value of the Knudsen parameter for which the WIMP energy transport 



is most efficient. We remark that other approaches are possible. In Refs. |L5| and |g9|, |30 



for example, a suitable interpolation of Eqs. flj]) and @ is used as a general expression for 
WIMP luminosity. We have checked numerically that our results are essentially independent 
of the chosen approach. 

As mentioned before, WIMPs can interact with nuclei either via coherent interactions or 
via spin-dependent ones. In the first case the WIMP cross section on the i-th nucleus of 
mass rrii is related to the corresponding WIMP-nucleon cross-section, a c , by: 

x 4 / i \ 2 



a t = a c (— pL_^i , (14) 
\m n J \m x + rrii J 

where m n is the nucleon mass. With this assumption, Eq. (||) may be rewritten as: 

1 - 1 E x jW r^) 3 ^4^)\ as) 



l x (r) l n (r)i \m n J \m x + m 

with: 

m n f l.Q ■ 10 -37 cm 2 \ / 150 g cm 



-3 s 



'-w = ^- — — • l-So- (16) 
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In the Sun the abundances Xi for nuclei heavier than Hydrogen and helium are much smaller 
than those of H and He. However, in the sum over elements of Eq.(|I5D, for heavy nuclei 
(mainly for iron) the reduction due to X$ is largely compensated by the coherence effect 
displayed by Eq.flnj). 

On the other hand, as far as spin-dependent interactions are concerned, the contribution 
of Hydrogen is overwhelmingly dominant over the contributions of other nuclear species 
with non-vanishing spins. Therefore, from now on we identify a s d with the spin-dependent 
WIMP-proton cross-section. 

In the following all WIMP matter interactions will be expressed in terms of a c and/or 
a s d- The notation a p will be used to denote a generic WIMP-proton cross-section; a p has 
to be identified with cr c (or a s d), whenever coherent (or spin-dependent) interactions are 
considered. 



III. WIMPS AND THE CENTRAL SOLAR STRUCTURE 

The efficiency of WIMP energy transport depends on several parameters. Specifically, it 
depends on the total number N x of WIMPs in the Sun, on the WIMP mass m x and on the 
WIMP-nucleon cross section a v . In order to determine the region of the WIMP parameter 
space where the WIMP energy transport is most efficient it is useful to define the quantity: 

* = lim^fc), (17) 
r->o L(r) 

where L is the radiative luminosity of the Sun. It is easy to show that: 

kJO) 



[Kj K ) 2 + 1 



«x(0) 



(18) 



where k 7iX (0) are the opacities at the solar center for radiation and WIMPs. 

In Figs. |l|-^| we show the lines in the plane (a p , N x ) which correspond to fixed values 
of S. We have chosen m x = 50 GeV as a representative value for the WIMP mass. Fig. 
[I] refers to coherent WIMP scattering, while Fig. ^ is obtained for spin-dependent WIMP 
interactions. The lines have been calculated by using the physical and chemical parameters 
of the Sun as they are obtained by a standard evolutionary code [j33fl . 

The qualitative behavior of the iso-5 lines is easily understood. The efficiency of the 
WIMP energy transport is proportional to l/er p (i.e. to K) in the conductive regime and 
to ij p (i.e. to K~ l ) in the non-local regime. The transition between the two regimes occurs 
at K = l x /r x ~ 0.4, i.e. at a cross-section a t which corresponds, for coherent scattering, 
to the value a t ~ 10~ 37 cm 2 and, for spin-dependent interactions, to o t ~ 10~ 34 cm 2 . The 
difference in the values of the transition cross sections at is due to the fact that, for the spin- 
independent case, interactions on elements heavier than Hydrogen is important, due to the 
coherence effect in the scattering which is manifest in Eq. (|i"4|), while for spin-dependent 
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interactions no coherence is present and therefore only scattering on Hydrogen (which is 
largely more abundant than other elements with non vanishing spin) matters. 

To be more quantitative, we can also show the analytic expression of the iso-5 lines as a 
function of o p . 

N x (6, a p ) = N{5) [(a p /a t ) + KM)" 1 ] , (19) 
where N{§) is given by the condition: 

K J (0) = 5K x (0;l x = 0Ar x ), (20) 

which gives: 

N(A\ An ,T3 f^T 1 

(21) 



N(S) ,4a rad T 3 / 
7T 3 / 2 r x = 3 °- 4r x A x 



and lj ~ 5 • 10~ 3 cm is the photon mean free path at the center of the Sun. 

We remark that, since WIMP radial distribution is rapidly decreasing (see Eq.@), the 
parameter 5 essentially provides an upper limit for the ratio L x /L in the Sun. As a conse- 
quence, the condition that S is not too small should be used as a necessary (but not sufficient) 
condition for WIMPs to have relevant effects on the Sun. This means that: 

i) We can safely conclude that regions of the parameter space for which 5 < 10~ 2 are 
not interesting. This point can be easily understood. WIMPs modify the relation between 
the luminosity and the gradient of temperature in the Sun. This essentially corresponds to 
a redefinition of the Sun radiative opacity k t The effect is maximal at the center of the 
Sun where one has /(l + 5). The radiative opacity is affected by uncertainties at 

the level of a few percent [33|. This means that we cannot disentangle WIMP effects if 5 is 



smaller than 10~ 2 . 

ii) We cannot conclude, at this stage, that regions in which 5 is large (say e.g. 5 > 1 ) 
correspond to sizeable modifications of the solar structure. The parameter 5 is, in fact, a 
local parameter. A large 5 only implies that WIMPs give a large contribution to the energy 
transport at the center of the Sun. In order to understand whether WIMPs have relevant 
effects on global properties of the Sun, we also have to compare the radius r x of the WIMPs 
extension region with the relevant length scales in the Sun (e.g. the dimension of the neutrino 
production region, the temperature height scale). 



IV. THE NUMBER OF WIMPS IN THE SUN 

Up to this point we have treated m x , a p and N x as independent quantities. However, 
this is not the noted long ago in Ref. [0. N x depends on the capture rate of the 

halo WIMPs by the Sun during its lifetime; in turn, this capture rate is proportional to the 
product of the WIMP density in the halo p x times a p , i.e. to the product £a p . Moreover, 
captured WIMP may annihilate with captured antiWIMP. Therefore N x also depends on 
the quantity (cr a v )o-> where a a is the WIMP-antiWIMP pair annihilation cross-section and 
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v is the relative velocity of the annihilating pair, while brackets indicate average over the 
WIMP velocity distribution in the Sun core. Explicitly, one has [BT]]: 



N x = T c t tanh(t Q /r), (22) 

where t & ~ 4.5 Gyr is the age of the Sun, T c is the WIMP capture rate and r = 1 / y/T c C a 
is a time scale parameter which determines whether and when equilibrium between capture 
and annihilation occurs. The quantity C a is given by: 

a ~ v20GeV ) ' 

with Vq ~ 2.3 x 10 28 cm 3 for the Sun. In the limit in which (cr a v)o —>■ 0, one has r — > oo 
and Eq. fl2"2"|) becomes: 

N x = T c t Q = N^ ax . (24) 

jymax re p resen |; S the upper limit for the WIMP number in the Sun at present time. 

The expression of the capture rate T c is rather involved. For a Maxwellian WIMP velocity 



distribution in the galactic halo it may be cast in the form |[35|| : 

1/2 



- 

<Ti— V 
. m x J 



esc ((h) 



^(00)^ , (25) 



where the cross sections Uj are related to o v as discussed in the previous Section, Mi is the 
total mass of the z-th nucleus in the Sun, v esc ~ 618 km sec -1 is the escape velocity at 
the Sun's surface and v = 270 Km s -1 is the WIMP velocity dispersion. The quantities 
(4>)i denote the reduced gravitational potential <p(r) = vl sc {r) / l v 2 esc averaged over the mass 
distribution of the i-th nucleus, as defined in Eq. ( |X^ ) in the Appendix. The quantity 
£(oo) ~ 0.75 is a suppression factor due to the motion of the Sun through the halo and the 
factor Si includes the effects of WIMP-nucleus mass mismatch and finite target dimensions 
(note that the <7; are point-like cross-sections). The expression of Si is given in the Appendix 
for completeness. The chemical composition of the Sun and the values of the quantities (</>);, 
which have been used in our calculations, are given in Table |. We note that for relatively 
heavy WIMPs like the ones which are considered in this analysis (m x > 30 GeV) the capture 
in the Sun for spin-independent interaction is not dominated by scattering on Hydrogen, 
even though this is the most abundant element, but instead by interactions on heavier 
elements. For instance, for a 50 GeV WIMP, capture is dominated by He, O and Fe. The 
smaller abundance of heavier elements is largely compensated by the coherence effect in the 
cross section shown in Eq. ([14]). On the contrary, in the case of spin-dependent interactions 
the only relevant capture process is on Hydrogen, since in this case no coherence effect is 
present and therefore capture is dominated by interactions with the most abundant element 
which possesses spin, i.e. Hydrogen. 
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Eq. (p5[) holds when the cross section a p is small enough so that multiple scatterings can 
be neglected. However, when a p is so large that every WIMP crossing the Sun is captured, 
T c saturates to a maximal value, and Eq. (p5|) is replaced by: 



where £ ~ 1.77. In Eq. (ffiSQ, the first term in the last square brackets refers to WIMPs 
whose orbit cross the Sun even without gravitational deflection, while the second term is a 
gravitational focusing factor and corresponds to those WIMPs whose orbits pass through 
the Sun because of gravitational deflection. One can easily see that for the Sun the focusing 
factor is dominant over the purely "geometrical" one, since the velocity dispersion of WIMPs 
in the Galaxy is a factor of two smaller than the escape velocity at the Sun's surface. 

The behaviour of the capture rate, and consequently of the total number N x of WIMPs 
captured by the Sun, as a function of the WIMP parameters is now easily understood. At 
fixed m x , N x increases linearly with o p up to the saturation level. For m x = 50 GeV, 
this behaviour can be seen in Fig. [l| (coherent interactions) and Fig. |] (spin-dependent 
interactions), where the thick solid line represents N x as a function of a v for the case of 
no-annihilation [o a = 0) and the thin solid lines show N x for some representative non- 
vanishing values of (cr a v) . In calculating N x as a function of <r p , we have set p x = pi, i.e. 
£ = 1. For pi we have used, here and throughout the paper, the default value pi = 0.3 
GeV cm -3 . The current uncertainty on p\ (0.2 GeV cm~ 3 < pi < 0.7 GeV cm -3 , for an 
isothermal galactic halo 0]) may actually increase or reduce N x by at most a factor of 
2. The no-annihilation line represents, for a given WIMP mass, the maximal number of 
WIMPs that can be captured by the Sun during its lifetime. If annihilation is present, N x is 
obviously smaller, depending on the strength of the annihilation cross section. Results very 
similar to those shown in Figs. [2] are obtained for different WIMP masses. Notice that the 
transition from the linear regime to the saturation level occurs at quite different values for 
the relevant WIMP-nucleon cross section in the two limiting cases of pure coherent and pure 
spin-dependent interactions. This is again due to the fact that for WIMPs heavier than tens 
of GeV, scattering on heavy nuclei in the Sun is enhanced for coherent interactions, and this 
shifts the transition from the linear to the saturation regime toward WIMP-nucleon cross 
sections smaller than in the spin-dependent case. 

The previous discussion may be easily extended to the case in which the WIMP under 
consideration does not provide the total amount of local dark matter (i.e. £ < 1). In this 
case, the relevant quantities are conveniently plotted in the plane ^a p -N x . For the case of 
coherent interactions the results are displayed in Fig. |3], whose four panels refer to different 
representative values: £ = 1, 0.1, 0.01, 0.001. In each panel the dashed curve denotes a 
representative iso-5 line with 5 = 1 and the solid line displays N x as a function of £cr p 
for the limiting case a a = 0. For the sake of comparison, the top-left panel repeats the 
case already shown in Fig. [1] (£ = 1). The features of the panels with £ < 1 are easily 




c + 



2v 2 



£(oo)7r#| 



(26) 
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understood in terms of those represented in the top-left panel: i) the iso— S curve simply 
shifts horizontally to the left by an amount equal to £, since 5 depends on the scattering cross 
section a v but it does not depend on £; ii) the slanted part of the solid line is unchanged 
since T c in the linear regime, as given by Eq. (|25|) , is a function of the product £<t p ; iii) 
the flat part of the solid line (saturation regime) is lowered by an amount £ since T c ^ sat 
in Eq.(P5|) depends on £ but is independent of a p . We explicitely consider the possibility 
of subdominant WIMPs since this is a feature which can naturally occur for specific dark 
matter candidates, like, for instance, the neutralino, as it will be discussed in Sect. [VI I[ 

We are now able to understand which regions of WIMP parameter space potentially lead 
to observable modification of the solar structure, by comparing the iso-5 lines with the 
curves which describe the number of WIMPs in the Sun. The first result is that the WIMP 
maximum effect occurs when annihilation is negligible and for scattering cross sections in 
the interval 10 -38 cm 2 < a c < 10~ 36 cm 2 for coherent interactions and 10 -35 cm 2 < <j s d < 
10~ 33 cm 2 for spin-dependent interactions. These cross sections correspond to a situation in 
which the number of WIMPs in the Sun has essentially reached its saturation value and the 
WIMP energy transport is maximally efficient (i.e. K ~ Kq). This is manifest in Figs. 0, |2| 
for m x = 50 GeV and occurs also for all the WIMP mass range m x > 30 GeV under study 
in the present paper. The value of 6 which corresponds to this maximal effect is 5 ~ 10, 
both for coherent and spin-dependent interactions. We will use this information in the next 
Sections where we will calculate the maximal deviation of the Sun properties when WIMPs 
are incorporated in its interior. 

We stress that, for the WIMP mass range considered in the present paper, the process of 
WIMP capture cannot produce modification of the Sun which refer to values of S larger than 
about 10, for any value of the scattering cross-section. Independent limits on the WIMP- 
nucleon interactions, like the ones that are obtained from direct search experiments, can 
further bound the values of 5 that can be reached. This is the case for coherent interactions, 
for which direct search experiments have already been able to set stringent bounds. For 
m x = 50 GeV, the upper limit on £cx c is shown in Fig. |l] by a vertical dashed line. From this 
figure, we can therefore see that once the direct detection constraint is taken into account, 
the largest values of 5 which can be obtained are of the order of 10~ 6 -10 -7 , which are far 
too small to lead to observable effects, as already discussed at the end of Sect. fTJ. Once 
rescaling is considered, the behaviours of the iso-5 and iV™ ax lines plotted in Fig. |3] show that 
<5's of the order of 10~ 3 -10~ 2 can be reached. These values are close to the necessary (but 
not sufficient) condition for WIMPs to have observable effects on the Sun (i.e. 5 > 10 -2 ), 
introduced in Sect. |T|, and therefore they need to be analyzed in detail. 

In the case of spin-dependent interactions, direct search experiments are also setting 
bounds on £,<J s d, but these limits are not model independent since they rely on the way the 
interaction of WIMPs with the spin of the nucleus is realized || . For the representative case 
of a Z-boson mediated interaction and for m x = 50 GeV, the present upper limit is shown in 
Fig. by a dotted vertical line. Since this line cannot be assumed to be a model independent 
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bound, in the case of spin-dependent interactions we are left with the possibility to reach 
values of 5 ~ 10. 

Up to this point we have quantified the amount of possible deviations induced by a generic 
WIMP in terms of the parameter 5. We now turn to discuss the possibility to experimentally 
observe such modifications by means of helioseismology and of the Boron solar neutrino flux. 



V. WIMPS AND HELIOSEISMOLOGY 

Helioseismology has provided very accurate information on the solar structure and it has 
been able to establish severe constraints and tests of the solar standard model (SSM) cal- 
culations. For instance, helioseismology accurately determines the depth of the convective 
envelope Rb and the photospheric helium abundance Y p h- Moreover, by inversion of helio- 
seismic data, one can determine the (squared isothermal) sound speed in the solar interior, 
u = P/p, with high accuracy (P and p denote the pressure and density inside the Sun). In 
Fig. H, we show the uncertainty in the helioseismic determination of u as a function of the 
radial coordinate R/R & . The light band corresponds to the so-called "3<r" errors which have 
been estimated conservatively by adding linearly all known individual uncertainties HT8"| . If 
uncertainties are added in quadrature, the global error is about one third. This yields the 
so called "la" errors which are shown by the dark band in Fig.|j. This latter procedure was 
also used by Bahcall et al. [[3(J with similar results. For intermediate values of the radial 



coordinate, the uncertainty in the squared sound speed is at the level of 10~ 3 , while in the 
innermost part of the Sun, as well as on the surface, it is at the percent level. We remark 
that, as found in Refs. [p| and [Q, uncertainties corresponding to statistical errors of the 
frequency measurements are generally much smaller than the "systematical uncertainties" 
of the inversion method. These latter include the choice of the reference model and of the 
free parameters in the inversion procedure. These uncertainties are particularly important 
close to the solar center, since there are very few p-modes that sample this region well. 
The errors quoted above include both statistical and systematic uncertainties. On the other 
hand, Refs. P7I and [38| present sound speed profiles with errors which look significantly 



smaller than our estimate. This occurs since "the quoted uncertainties include only the 



contribution arising from the frequency observations", as clearly stated in [38]. 

Fig. f| also shows the difference between u as predicted by the SSM of Ref. |39| and 
the helioseismological determination u sun , normalized to u sun . We notice that the SSM 
accurately reproduces u for all the radial profile. 

In order to understand whether helioseismology can tell us something about WIMPs, we 
have constructed solar models for values of the WIMP parameters which maximize the role 
of WIMPs in the energetics of the Sun, irrespective of possible bounds derivable from direct 
search experiments. As discussed in the previous Section, this corresponds to 5 ~ 10. For 
m x = 50 GeV, this is shown, for instance, by the large dot on the "no-annihilation line" in 
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Fig.[l]. We have taken into account the time evolution of the number of WIMPs in the Sun 
and we have described WIMP energy transport at each stage of Sun evolution according to 
Eq. fll2j) . The sound speed profiles obtained for the representative WIMP masses m x = 30, 
50 and 100 GeV are compared with the sound speed profile of the SSM in Fig.|5]. Larger 
WIMP masses lead to even smaller effects. For any WIMP mass, smaller values of 5 clearly 
correspond to smaller effects. 

The first observation is that, even in the most favourable situation, WIMPs produce a 
modification of the sound speed profile which is much smaller than the accuracy of present 
helioseismic determinations. The smallness of the effect is explained by two main reasons. 
First, WIMPs with masses larger than tens of GeV are confined in a small region of the Sun of 
radius R < O.O2i? . In this region the relevant physical quantities (temperature, luminosity, 
etc.) have small variations. As a consequence, even if WIMPs affect sizeably the gradient of 
these quantities (e.g. the temperature gradient), this translates into small modifications of 
their "global" radial profile. Second, a compensation occurs between two different effects. 
This can be easily understood in terms of the fact that, from the perfect gas law, which 
describes to a good approximation most of the solar core, one has u = P/p ~ T/fi, where 
fi is the mean molecular weight. WIMPs produce a quasi- isothermal core, thus decreasing 
the central temperature T c of the Sun. At the same time, however, with a lower T c one 
has smaller nuclear reaction rates. This translates into a slower Helium production at the 
center of the Sun and, thus, into a decrease of the mean molecular weight. When considering 
u ~ T/n, the effects of the variations of T and /i compensate, leaving the sound speed profile 
almost unaffected. 

The results shown in Fig. |5| refer to a situation which maximizes the effect of the presence 
of WIMPs in the Sun: 5 ~ 10. As discussed in the previous Section, for spin-dependent 
interactions, this value of 5 can be reached, while for coherent interaction <5's of the order of 
at most 10" 3 -10- 2 can be obtained. Our results show that in neither case (spin-depenent 
nor coherent) the presence of WIMPs in the Sun can affect u to a level which is currently 
accessible by helioseismology. 

We notice also that the necessary (but not sufficient) condition stated in Sect. [TTJ, i.e. 
that 5 must be larger that about 10~ 2 in order to possibly have an observable effect, may 
be reinforced: in fact our detailed calculation shows that, for WIMP masses above 30 GeV, 
not even 5 ~ 10 is sufficient to produce a currently observable effect. 

In summary, the results of this Section show that no significant information about WIMPs 
with masses m x > 30 GeV can be obtained at present from helioseismic data, for any values 
of the WIMP-nucleon cross section. 

Conclusions conflicting with ours were derived in Ref. p5| . This may be traced back 



to the fact that in Ref. |29|] the accuracy of the squared sound speed is taken to be of 



about 0.1% also in the innermost part of the Sun core relevant for WIMP confinement, say 
R < 0.02R Q . We argued above that, at these small radii, the accuracy of the helioseismic 
determination is significantly worse. 
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VI. WIMPS AND THE BORON NEUTRINO FLUX 



From the previous discussion, one could expect that the Sun central temperature T c is 
more sensitive than the central sound speed to the presence of WIMPs in the Sun. This is 
actually the case, as can be seen in Table ||, which shows that WIMPs produce variations 
of T c at the percent level, while the central sound speed u c is affected at most at the level 
of 1(T 3 . 

Unfortunately, the central temperature of the Sun is not directly observable. The most 
direct information on it comes from the the measurements of the solar neutrino fluxes. In 
particular, the comparison between the SNO charged current and Super-Kamiokande data 
|2"T| , |2"2] , |2"B| . |2"4] , |2"5] , P^ ] has allowed to determine, in a model independent way, the 8 B neutrino 
flux produced in the Sun with an accuracy of about 20% 2 . This result, under the assumption 
that the 8 B neutrino flux scales with the Sun central temperature as <3>b oc T c 20 , can be used 



to determine T c at the 1% level [[28 



However, the presence of WIMPs in the core of the Sun produces a very peculiar modifica- 
tion of the temperature profile for which the scaling law $b oc T c 20 is not valid. WIMPs can 
sizeably affect the temperature profile only in the innermost part of the Sun, and in partic- 
ular in a region which is smaller than the Boron neutrino production region (tb — 0.05R & ). 
Fig.|5] shows the radial profile of the Sun temperature in the SSM and in models of the Sun 
where WIMPs with masses of 30, 50 and 100 GeV are present. Again, we have chosen the 
WIMP parameters which maximize the effect of their presence in the Sun. Fig.|6] clearly 
shows that even though the central temperature of the Sun can be lowered by the presence 
of WIMPs at the percent level (as it is reported also in Table [Tl[), nevertheless the temper- 
ature at rs — 0.05-R©, where the Boron neutrino production is maximal, is left practically 
unchanged. 

In order to perform a quantitative analysis, one has therefore to look directly at the 
Boron neutrino flux and not simply at the central temperature T c . If this is done, it turns 
out that the solar models of Fig. ||, which correspond to WIMP masses equal to 30, 50 
and 100 GeV, predict variations of the Boron neutrino flux equal to: 5§b/&b — —3.6%, 
— 1.0% and —0.1%, respectively. These variations are well below the accuracy of the present 
determination of the Boron neutrino flux. We remind that, for each WIMP mass, we have 
optimized the other WIMP parameters in order to maximize the effect. Larger masses and 
WIMP annihilation would lead to even smaller deviations. 

In conclusion, our results allow us to conclude that no information about WIMPs with 



2 The SNO experiment has recently reported the observations of neutral current v interactions on deuterium 
[f27| confirming the 8 B neutrino flux predicted by SSM with an accuracy of about 12%. This result is 
obtained under the assumption that the solar neutrino energy spectrum is undistorted. For this reason, in 
our analysis, we refer to the result obtained by comparing SNO charged current and Super-Kamiokande 
data, which instead do not rely on this assumption, see Refs. ElL M. 
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massess m x > 30 GeV can be obtained from our present knowledge of the Boron neutrino 
flux. 

A previous independent analysis |30| concluded that WIMPs lighter than 60 GeV are in 
conflict with the present accuracy in the determination of the Boron neutrino flux, based 
on the fact that the presence of WIMPs in the core of the Sun can modify the central 
temperature T c at the 1% level. We should stress that this lower limit on the WIMP 
mass does not apply, since, as we have discussed above, WIMPs, even though capable of 
changing T c at the percent level, nevertheless are too much concentrated in the interior of 
the Sun to produce sizeable modifications of the Boron neutrino flux, which is the quantity 
experimentally measured. 



VII. ONE REALIZATION OF WIMP: THE NEUTRALINO 

It is instructive to discuss how some of the previous properties are realized in the case of 
one of the favourite WIMP candidates: the neutralino. In Figs. [7]-|9] we show some of the 
relevant properties. 

The supersymmetric scheme we are adopting here is an effective Minimal Supersymmet- 
ric extension of the Standard Model (MSSM). Its parameters are defined directly at the 
electroweak scale and represent the minimum set necessary to shape the essentials of the 
theoretical structure of the model and of its particle content. We refer to Ref. ET| for details 
on the theoretical aspects and on the updated experimental bounds. 

The neutralino-nucleon cross sections, both coherent and spin-dependent, are calculated 
in terms of the couplings and supersymmetric particle masses as discussed, for instance, in 
Refs. P, ITU]. The parameter £ is evaluated by using a standard rescaling procedure Wl 



which relates the fraction of local neutralino dark matter in the Galaxy to its fractional 
abundance in the Universe, i.e. to its relic abundance fl x h 2 : 



£ = min 



n x h 2 



(n m h 2 ) 



(27) 



where (f2 m /i 2 ) m m is the minimum value of Q m h 2 compatible with halo properties. Here 
we set (fi m /i 2 )min = 0.05. The present-day neutralino relic abundance has been calculated 



according to Ref. [43 



To establish some relevant order of magnitude for the annihilation cross-section o a it 
is useful to recall the approximate expression: Q x h 2 ~ 3- 10 -27 cm 3 s~ 1 /(cr a v)i nt , where 
(c>a v)i n t denotes the thermal average of the product (cr a ■ v) integrated from the freeze-out 
temperature Tj to the present-day one. For instance the upper bound Q x h 2 < 0.3 translates 
into a lower bound (<r a v) int > 10~ 26 cm 3 sec -1 . 

In Fig. ^ we display the correlation between £ and (<r a v)o, calculated by varying the 
MSSM parameters and applying to the model the relevant experimental constraints. Fig. [7] 
shows that for a neutralino in the MSSM, the rescaling factor is always smaller than 1 for 
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values of (<r a v)q larger than about 10 _27 -10~ 26 cm 3 s -1 . Notice that (a a v)o denotes the 
thermal average of the product (a a ■ v) in the Sun's core, and does not coincide with the 
relevant thermal average which enters in the calculation of the relic abundance. 

We turn now to the graphic representations in the plane £,a p -N x . Fig. || (Fig. |fj) shows 
the scatter plot of the MSSM configurations for which the total capture rate is dominated 
by coherent (spin-dependent) interactions and for the neutralino mass range: 50 GeV < 
m x < 500 GeV, the lower bound arising from experimental constraints on the MSSM. The 
points are labelled according to the values of the neutralino relic abundance: dots denote 
configurations where the neutralino is the dominant dark matter component: 0.05 < £l x h 2 < 
0.3, i.e. £ = 1; crosses refer to configurations where the neutralino is a sub-dominant dark 
matter component: Q x h 2 < 0.05, i.e. £ < 1. The slanted solid line shows N x for m x = 50 
GeV and for the limiting case of non-annihilating WIMPs. This line represents, for each of 
the two limiting cases of coherent and spin-dependent domination, the actual upper limit 
for N x for neutralinos in the effective MSSM, due to the lower bound on the neutralino mass 
quoted above. The dot-dashed lines denote iso-<5 contours for two representative values of 
5. For each pair, the upper line refers to m x = 50 GeV and £ = 1, while the lower line refers 
to m x = 500 GeV and £ = 1. 

By comparing the scatter plot with the iso-5 lines, we can notice that the configurations 
with no rescaling (£ = 1) are compatible with values of 5 smaller than about 10~ 9 , both 
for coherent and spin-dependent interactions. We remind that these are configurations 
with large neutralino relic abundance. In the case of configurations with £ < 1, we have 
to consider the shift of the iso-5 lines shown in Fig. |3| and the range of values for £ in 
the MSSM shown in Fig. [7[ In this case, we obtain that the maximal effect induced by 
subdominant relic neutralino is for values of S of the order of 10 -6 for the coherent case and 
10~ 4 for the spin-dependent one. 

We can therefore conclude that the maximal effects induced in the Sun by relic neutralinos 
are much smaller than those valid for a generic WIMP and discussed in Sect. |TV|. This implies 
that the specific case of the neutralino has even smaller chances to be bounded by the study 
of solar properties. 

VIII. CONCLUSIONS 

Our understanding of solar physics has significantly advanced in recent times, because 
of remarkable improvements in helioseismology and of new data on solar neutrinos. In the 
present paper we have addressed the question of whether these new developments put solar 
physics in a position to provide constraints on the possible presence of WIMPs in the core 
of the Sun, for values of WIMP parameters under current exploration in WIMP searches. 
We summarize here our main results: 

• We have provided a quantitative criterium to determine whether a putative WIMP 
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candidate could produce observable modifications of the solar structure. Namely, we 
have introduced the parameter 5, defined as the ratio between the WIMP luminosity 
and the radiative luminosity at the center of the Sun (see Eg. (|T7|)), which can be 
analytically calculated as a function of the WIMP-nucleon scattering cross section o v 
and the number of WIMPs in the Sun N x . By considering the uncertainty in the 
radiative opacity, one finds that 5 > 10~ 2 is a necessary (but not sufficient) condition 
for WIMPs to have observable effects on the Sun. 

• We have calculated the number of WIMPs in the Sun for a generic WIMP candidate 
as a function of the WIMP-nucleon scattering cross section a p and of the WIMP- 
antiWIMP pair annihilation cross-section a a , both for coherent and spin-dependent 
WIMP interactions. We have shown that a value 5 ~ 10 can be reached, if WIMP- 
antiWIMP annihilation is negligible and a p ~ 10 -37 cm 2 (coherent scattering) or 
a p ~ 10~ 34 cm 2 (spin-dependent interactions). For coherent WIMP scattering, such 
large values for the cross sections are excluded by direct search experiments, unless 
WIMPs give a subdominant contribution to dark matter of our galaxy. In this case, if 
one takes into account direct search bounds (£cr p < 1CT 41 cm 2 , where £ is the rescaling 
parameter which accounts for the actual fraction of local dark matter to be ascribed 
to the given WIMP), one obtains 5 < 10~ 6 in the assumption of no rescaling (£ = 1) 
and 5 < 10~ 3 in the general case (£ < 1). 

• We have considered the neutralino as a specific WIMP candidate. We have calculated 
the number of neutralinos in the Sun for MSSM configurations compatible with present 
experimental and cosmological constraints. For each configuration, we determined the 
rescaling parameter £ by the standard rescaling procedure described in Eq. (j27|). In the 
case of no rescaling (£ = 1), we obtained 5 < 10~ 9 . For subdominant neutralinos (£ < 
1), we obtained 5 < 1CT 6 for configurations in which the capture rate is dominated by 
coherent scattering and 5 < 10~ 4 for configurations in which the dominant contribution 
is due to spin dependent interactions. 

• The previous points already show that solar physics is not competitive with direct 
experiments in a large part of the WIMP parameter space under current exploration 
in WIMP searches. In order to complete our analysis and to understand whether some 
information could be obtained from the present helioseismic and solar neutrino data, 
we have constructed solar models, for WIMP masses above 30 GeV, choosing the value 
of the WIMP parameters which maximize the effect of WIMPs on the Sun (5 ~ 10). 
As a result of the presence of WIMPs in the Sun, we obtained variations of the sound 
speed profile and of the 8 B neutrino flux which are within the current experimental 
uncertainties. The smallness of the effects is essentially due to the the smallness of 
the WIMP estension region which, for WIMPs with masses larger than tens of GeV, 
is r v < 0.02i? M . 
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In conclusion, no constraints can be derived at present from solar physics for WIMPs 
with masses above 30 GeV. Our conclusions are at variance with results derived in Refs. 
( p9| , p0|1 ). The origins of these disagreements have been elucidated in the present paper. 



APPENDIX A 



We give here for completeness the expression of the factor Si, introduced in Eq . (P5|) , as 
derived from Ref . [|35|j . It is given by: 



S; 



{4>)i CM 



(Al) 



where the brackets indicate an average over the mass density profile Pi(r) of the i-th nucleus 
in the Sun: 

{f)i = ^rf ^r 2 dr Pi (r) f(r) (A2) 
Mi Jo 

All the quantities are defined in Section III with the exception of Fi, which is given by: 



vie 3b v 



X(-V,V) - x(A-,A+) 



exp[—ar]) 
(l + a)V2 



XK V) - X(A-, A + )} 0£ exp[-(a - b)A*} 



(A3) 



The definition of the quantities a, b, r], A, A+_, A+-, fj, and fj may be found in Ref. |35|| . 
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TABLE I: Sun composition, average mass fractions of the chemical elements in the Sun as obtained 



from Ref.|44| and values of the average reduced gravitational potential {<j))i. 
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TABLE II: Fractional variations of the central temperature of the Sun T c , of the squared isothermal 
sound speed in the Sun interior and of the 8 B neutrino flux, calculated in solar models where WIMPs 
are accreted in the Solar core. For each model, the WIMP-nucleon scalar cross section has been 
set at the value which maximizes the effect of the presence of WIMPs in the Sun. 



m x (GeV) 


ST C /T C 


5u c /u c 




30 


-2.1 • 10" 2 


-4 • 10~ 3 


-4.5 • 10" 2 


50 


-1.2 • 10~ 2 


-2 • 10 -3 


-1.3- 10~ 2 


100 


-5.2 • 10~ 3 


-6 • 10~ 4 


-2.1 • 10- 3 
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A. Bottino, G. Fiorentini, N. Fornengo, B. Ricci, S. Scopel, F.Villante (2002) 




FIG. 1: The solid (red) curves denote the total number N x of WIMPs accumulated in the Sun's 
core due to gravitational capture as a function of the WIMP-nucleon cross section, for a WIMP 
mass of 50 GeV and for the case of purely coherent WIMP-nucleus interactions. The thick (red) 
solid line refers to WIMPs which do not annihilate once captured. The thin (red) solid lines 
refer to annihilating WIMPs, for increasing values of the zero-temperature thermally-averaged 
WIMP-antiWIMP annihilation cross section: {cr a v)o = 1CT 30 , 1(T 27 , 1CT 24 , 1CT 21 cm 3 s -1 . The 
dot-dashed (blue) lines show the contours of constant values of the parameter 5, defined in Eq. fll7|) , 
which quantifies the efficiency of energy transport of WIMPs with respect to photons at the center 
of the Sun. The vertical (green) dashed line denotes the upper limit on the WIMP-nucleon elastic 
cross section for coherent interactions. The large (red) dot identifies the WIMP parameters for 
which the WIMP effect on the Sun's energetic is maximal. 
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A. Bottino, G. Fiorentini, N. Fornengo, B. Ricci, S. Scopel, F.Villante (3002) 




FIG. 2: The same as in Fig. [l], for spin-dependent WIMP-nucleus interactions. The vertical (green) 
dotted line denotes the upper limit on the WIMP-nucleon cross section for Z-boson mediated spin- 
dependent interactions. 
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FIG. 3: The same as in Fig. H, plotted as a function of £a p , for different values of the rescaling 
parameter £. The (red) solid lines denote the total number of WIMPs accumulated in the Sun 
for a non-annihilating WIMP. The (blue) dashed lines show the contours of iso-<5 for 6 = 1. The 
vertical (green) dashed lines indicate the upper limit on ^a p for coherent interactions. 
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FIG. 4: Uncertainty bands in the reconstruction of the squared isothermal sound speed u in the 
Sun interior. The light (yellow) region corresponds to the conventional "3cr" error band, while the 
dark (blue) band identifies the conventional "la" range. The solid line corresponds to the fractional 
difference between u determinaed from observational data and u calculated from the BP2000 solar 
model H). 
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FIG. 5: Relative difference of the squared isothermal sound speed in the Sun interior u calculated 
for a standard solar model and for solar models with accreting WIMPs, as a function of the radial 
coordinate R/Rq. The light (yellow) region corresponds to the conventional "3cr" error band, while 
the dark (blue) band identifies the conventional "lcr" range. The dotted, dot-dashed and dashed 
lines refer to solar models with accreting WIMPs of masses m x = 30, 50, 100 GeV and coherent 
scattering cross sections which, for each mass, maximize the effect of the presence of WIMPs in 
the Sun. 
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FIG. 6: Temperature profile of the standard solar model and of models of the Sun where WIMPs 
with masses of 30, 50 and 100 GeV are present. The coherent WIMP-nucleus scattering cross 
sections are chosen, for each mass, in order to maximize the effect of the presence of WIMPs in 
the Sun. 
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FIG. 7: Rescaling parameter £ = min(l, Q x h 2 /0.05) as a function of the zero-temperature 
thermally-averaged neutralino self-annihilation cross section {cr a v)o in the effective minimal su- 
persymmetric standard model (MSSM). 
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Capture: coherent dominated 

A. Hotlino, G. Florentine, N. Pornengo, B. Ricci, 3. Scopel, F.Villante (2002) 
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FIG. 8: Total number N x of neutralinos accumulated in the Sun as a function of £<7 C in the effective 
MSSM. Only configurations where the total capture rate is dominated by coherent interactions are 
displayed. Neutralino masses are in the range 50 GeV < m x < 500 GeV. (Red) dots denote 
configurations where the neutralino is the dominant dark matter component: 0.05 < ^ x ^ 2 < 0.3 
(i.e.: £ = 1). (Blue) crosses refer to configurations where the neutralino is a sub-dominant dark 
matter component: S7 x /i 2 < 0.05 (i.e.: £ < 1). The slanted (pink) solid line shows N x for m x = 50 
GeV and for the limiting case of non-annihilating WIMPs. The (black) dot-dashed lines denote 
iso-<5 contours. For each pair, the upper line refers to the iso-5 curve for m x = 50 GeV and £ = 1, 
while the lower line refers to the same iso-5 and £, for m x = 500 GeV. The vertical (green) dashed 
lines indicate the upper limit on £cr p for coherent interactions. 
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Capture: spin dominated 



A. Hotlino, G. Florentine, N. Pornengo, B. Ricci, 3. Scopel, F.Villante {2002) 




Log 10 (fa sd /cm 2 ) 

FIG. 9: The same as in Fig. B for MSSM configurations where the total capture rate is dominated 
by spin-dependent interactions. 
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